Quantum bright soliton in a disorder potential 
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At very low temperature, a quasi-one-dimensional ensemble of atoms with attractive interactions 
tend to form a bright soliton. When exposed to a sufficiently weak external potential, the shape of 
the soliton is not modified, but its external motion is affected. We develop in detail the Bogoliubov 
approach for the problem, treating, in a non-perturbative way, the motion of the center of mass of 
the soliton. Quantization of this motion allows us to discuss its long time properties. In particular, 
in the presence of a disordered potential, the quantum motion of the center of mass of a bright 
soliton may exhibit Anderson localization, on a localization length which may be much larger than 
the soliton size and could be observed experimentally. 

PACS numbers: 03.75.Lm,72.15.Rn,05.30.Jp 
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I. INTRODUCTION 

Anderson localization is a localization effect predicted 
to take place for a wave propagating in a disordered po- 
tential [1]. It is due to multiply scattered waves from 
random defects and yields exponentially localized den- 
sity profiles, resulting in a complete suppression of the 
usual diffusive transport associated with incoherent wave 
scattering [2]. While in the three dimensional world, one 
may observe a transition between extended and localized 
states, in a one-dimensional (ID) world, Anderson local- 
ization is a typical feature of the motion in a disordered 
potential [3|. 

Cold atoms form a wonderful toolbox for controlling 
parameters of the system under study [ij . It comes out 
as no surprise that attempts have been made for a direct 
observation of the Anderson localization in cold atoms 
settings. Already the first attempts ||, @, 0, i, H 03 
have revealed that the presence of atomic interactions 
may deeply affect the physics of the problem and make 
the observation of the localization non trivial. Further 
theoretical studies [HI, [H, [Hj] were followed by successful 
observations of the phenomenon made possible by going 
to the regime of very weakly interacting particles [3]. 
While in that work a random speckle potential was used, 
in another attempt fl5j a quasi-periodic version of the 
potential using superposition of laser beams was created 
resulting in the observation of Aubry- Andre [lril ] localiza- 
tion for noninteracting atoms. 

Anderson localization is a one-body phenomenon, and 
it is important to understand how it is modified when 
interactions between particles - in our case, cold atoms - 
are taken into account. In the absence of any external po- 
tential, at zero temperature, ID particles interacting at- 
tractively tend to cluster together, forming a bright soli- 
ton. Explicit solutions of the many- body problem can be 
found for a contact interaction [17] . Altogether, a bright 
soliton appears as a composite particle, whose position 
is given by the center of mass of the constituting atoms 
and a mass equal to the sum of the mass of the atoms 



(see next section). Using external potentials, it has been 
experimentally shown how to put solitons in motion [l8| . 
The purpose of this contribution is to discuss what hap- 
pens to a bright soliton exposed to a weak and smooth 
disordered potential [lj| [2(|. Of course, if that potential 
was sufficiently strong, it could probably destroy the soli- 
ton altogether, break it into pieces etc. We are, however, 
interested in the other limit when the external potential 
is sufficiently weak and smooth not to perturb the soli- 
ton shape. It is then quite reasonable to expect that, if 
this weak potential is of random nature (disorder) that 
the soliton as a composite particle, undergoes multiple 
scattering, diffusive motion and eventually Anderson lo- 
calization. In a recent short contribution [2l| we have 
shown that this is indeed the case by considering the ef- 
fective quantum motion of the soliton. The present work 
brings a detailed derivation of the effective Hamiltonian 
applied before, and shows examples of the corresponding 
localized eigenstates. It provides thus a complementary 
material to our previous work [2l| . 



II. MEAN FIELD DESCRIPTION 



Equations of motion for a bright soliton in a 
disorder potential 



Consider an ensemble of cold atoms (bosons) with 
attractive interactions at zero temperature. We as- 
sume a strong harmonic transverse confinement so a one- 
dimensional approximation can be used. In the mean 
field approach, a c-number function <j> takes the place of 
the bosonic field operator ip. is a solution of the Gross- 
Pitaevskii equation 
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where we have adopted the following natural units for 
energy, length and time, respectively 

Eq = 4mw^_o 2 , (2) 
2\a\mtu± 

h = — — -. (4) 



4a 2 moj^_ 

The transverse harmonic confinement frequency is de- 
noted by a is the atomic s-wave scattering length, 
and m the mass of an atom. We normalize <j) to the total 
number of particles N . Eq.{T]) admits a stationary bright 



soliton solution e~ 




where 



2£ cosh[(z - q)H\ 



(5) 



the chemical potential fi = —N 2 /8 and the soliton width 
£ = 2/N. This bright solitonic solution minimizes the 
energy functional 



E 



dz 



(6) 



Observe that eq. (O allows for an arbitrary center-of- 
mass (CM) position q and an arbitrary global phase 0. 

Suppose the soliton is placed in a weak and smooth dis- 
order potential, V(z), with variance Vq and correlation 
length (Jq. We will concentrate on the case when oo < £ 
but the approach we present is general. Linearization of 
the Gross-Pitaevskii equation allows us to describe the 
perturbation of the soliton due to the presence of a weak 
potential [23(. Indeed, the substitution 



(7) 



into ([TJ) supplemented with the potential V(z) leads to 
the following inhomogeneous time-dependent Bogoliubov 
equations 
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S = V(z) Uz-q). 



(9) 



(10) 



In Eq. ((8]) we have neglected terms of order higher than 
0(5(j), V). Solution of (|8|) can be expanded in right 
eigenvectors and corresponding adjoint modes of the non- 
hermitian operator C However, this operator is not diag- 
onalizable [23|, [24|, [2|| . For all eigenvectors (u n , v n ) corre- 
sponding to non-zero eigenvalues E n , the adjoint modes 
are left eigenvectors of the C. That is no longer true for 



the zero-eigenvalue modes. There are two zero modes in 
our system 



Ug 



= idff 



= id n 



(11) 



which are related to a small modification of the global 
phase of the solution ([5]) and to a small shift of the CM, 
respectively pll . [25j . As both modifications cost no en- 
ergy they appear as zero modes of the C operator. In- 
deed, it is consistent with quadratic expansion of the en- 
ergy functional, 



E = const + i J dz (5<f)* , - 



C 



(12) 



where we see that contributions to soliton perturbation 
from zero modes do not change E. The modes adjoint to 
the zero modes are 
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which has been found by solving 
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where Mg and M„ are determined by the requirements 
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1 and 
It turns out that 



uf\u q 



{vf\v q 



M n = N. 
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The latter is equal to the total mass of the system. Equa- 
tion (Q~||) ensures that (ug d , Vg d q ) are orthogonal to all 
eigenvectors of C with E n f^{). 

Perturbation of the soliton can be expanded in the 
complete basis vectors 



ug 
Vg 



Pa 



,,ad 
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q - q 



p„ 



, ad 
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n,E n >0 



,(16) 



where real q' and 0' describe translation of the soliton 
and shift of its global phase, respectively, while P q and Pg 
(also real) are momentum of the CM of the soliton and 
momentum conjugate to the global phase, respectively. 
The momentum Pg = N' — N represents deviation from 
the average total number of particles N. Deformation of 
the soliton shape is described by complex variables b n . 
Substituting (fTB)) into ([5]) and projecting on the basis 
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vectors results in a set of equations 
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d t q' = 
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-J dz \Mz-q)\ 2 9 z v(z), 
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wit eve real-valued 






s 


n = (u n \S) + (v n \S*). 


(22) 



Equation (|T7|) describes linear evolution of the global 
phase and it is possible to obtain 9'{t) = 6 = const by 
a proper choice of Pg. The latter is a constant of mo- 
tion, see (fl~8|) . We consider a weak disorder potential 
when do < £. Therefore the force acting on the CM, 
which is the force acting on a single particle convoluted 
with the soliton profile (f2"0"|) , is small and it oscillates 
around zero as a function of q. Thus, Eqs. (|19|) - (f20|) 



imply that, if we choose P q (0) = and such a q that 
/ dz\(j>o(z — q)\ 2 d z V = 0, then q'(t) — q — const. 



B. Deformation of the soliton shape 

We have seen that in a disorder potential the CM of the 
soliton can be fixed and its global phase can be constant. 
Let us now concentrate on the set of Eqs. (|2~lj) which de- 
scribe changes in the soliton shape due to the presence of 
a disorder potential. Solving Eqs. (|21[) with an assump- 
tion that initially the bright soliton is unperturbed, i.e. 
b n {0) = 0, we obtain 



Sr. 

En 

n,E n >0 
i I „iE„t 



= E £ [(«-**•' -i) 

l)v* n (z-q)]. (23) 



The lowest energy of the Bogohubov modes in the case of 
the bright soliton is E x = = N 2 /8 [26]. Thus a large 
gap in energy separates the soliton from the Bogohubov 
modes. These modes are delocalized and describe radi- 
ation of the soliton. The energy spectrum can be well 
approximated by a shifted free particle dispersion rela- 
tion 



2tt 2 



(24) 



where n is integer and L stands for the size of a box 
in which we consider our system. Moreover, due to the 
radiation character of the modes 



< 



< 




(25) 
(26) 



The latter inequality is obtained taking a rectangular 
profile of size £ for the bright soliton. Finally, with 
sin 2 (£'„t/2) < 1 and £ n l/En ~ / dn/E n , for defor- 
mation of the soliton shape, 



(27) 



(28) 



|0o + ^| 2 «|^o| 2 + 0o S^+^o S4>, 

we obtain the following estimate 

|0o 5<t>* + <j>t 5<j>\ <4|F |, 

and if it is much smaller than |</>o| 2 < 2|/x|, the shape 
of the soliton is negligibly changed. Hence, if we want 
the shape of the bright soliton to be unaffected by the 
presence of a disorder potential a sufficient condition is 



\v \ « H- 



(29) 



Note that the upper bound on Vq requires the potential 
to be sufficiently smooth, in particular the case of a 8- 
correlated disorder potential is excluded by this condition 

EH- 



C. Dziarmaga approach 



In Sec. Ill Al equations of motion for a bright soliton 
in the presence of a weak disorder potential have been 
obtained using the perturbative expansion (p~6|) . Conse- 
quently the long time evolution of the CM of the soliton 
for P q (0) 7^ cannot be described by these equations. 
Indeed, after a finite time \q'(t) — q\ > £ and the per- 
turbative approach breaks down. Similar problem may 
occur in the case of the 8' variable. 

We will be interested in a quantum description of the 
bright soliton where states corresponding to superposi- 
tion of the CM position over a distance much larger than 
£ will be considered. Therefore we need a method that al- 
lows us to describe non-perturbative displacement of the 
soliton. To this end we adopt Dziarmaga approach intro- 
duced in a problem of quantum diffusion of a dark soli- 
ton . Following Ref. [25| we do not perform a linear 
expansion of a perturbed soliton wave-function around 
fixed q and 9 like in (fl"6)) but we treat q and 9 themselves 
as dynamical variables 
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(30) 



Note that now if q(t) and 9{t) are changing in time all 
modes also evolve because they depend on q and 9, e.g. 
Ug d = u e td {z — q(t)). Substituting ([3H]) into energy func- 
tional © supplemented with the J dzV\(j)\ 2 term and 
keeping terms of order 0(P 2 ,b 2 , PV, bV) only, we obtain 
the effective Hamiltonian 



H 



P 1 

Q 

2M„ 



dzv(z) \M*-qf 



4 



p2 

2M 9 



■2P e {d N <h\V<fo) 



J2 (E n b* n b n + (b n + b* n )s n ), (31) 



n,B n >0 

which generates the following equations motion 

dH _ Pg 

dP e Me 
dH 
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dH 

dP~ q ~ 
dH 
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dH 

db* ~ 
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M a 



E„ b n 



(32) 
(33) 
(34) 



dz \cj) {z-q)\ 2 d z V(z), (35) 
+ s n . (36) 



In 

(bn 



we have neglected terms Pgd q {dN(f>o\V 4>q) an d 
b^dqSn because they are of order of 0(PV, bV) 
while in the equations we keep the linear terms only. 
Strictly speaking in order to show that the pairs of vari- 
ables in (f32"|) - ([36"|) are canonically conjugate one should 
switch to Lagrangian formalism of the problem, however, 
as the result is obvious, we have skipped it, see [25| . 

Equations (|3"2" ]) -(|3"S ]) possess a form identical to (fT7|) - 
(|2ip . However, q and 9 present in cf>o and s n on the 
right hand side of the current equations are not fixed 
and evolve in time. It introduces couplings between q 
and 9 and b n degrees of freedom which were absent in 
([i"7 ]) -([2T ]) . Inserting solutions of (32] ) -([36] ) into ([30]) we 
can obtain long distance propagation of a bright soliton 
including possible changes of its shape, something not 
possible with the expansion (TIT))) , 

The Hamiltonian (|3ip cannot be used for extremely 
large momentum of the CM. That is, it is valid provided 
P q £;/N < 1, compare (J3DD and ([15]). For the case of 
large P q see (ljjj . Note also, that due to the fact Mg 
is negative, see ([TS"]) . the bright soliton ([5]) is a saddle 
point of the energy functional ([5]). It has, however, no 
consequences since Pg = N' — N is a constant of motion. 



2Mg 



2Pe(d N <fi \Vct>o) 



J2 [EJlk + {b n + bl)s n ) , (37) 



where 



and 



n,E n >0 



P q = -id q , 



Pg = N-N = ~ldg, 



q,P q 


= i, 


9,Pg 


= i, 


, b m 





(38) 
(39) 

(40) 
(41) 
(42) 



Because [Pg,H] = we can choose a state \N) of 
the many body system with exactly N particles where 
Pg\N) — 0. If we consider the Bogoliubov vacuum 
state of the quasi-particle operators, i.e. & rl |0b) = 0, 
such a state will be very weakly coupled to other eigen- 
states of the E n b n b n operator because the coupling 
strengths s n are, for a weak disorder potential, much 
smaller than the large energy gap for quasi-particle ex- 
citations Ei = \fi\ = N 2 /8 [201 . Hence, the effective 
Hamiltonian that describes the CM motion reduces to 



(N;0 b \H\N;0 b ) 



P 2 

q 

2N 



dzV(z) |<M*-9)| 5 



(43) 



where we have inserted explicit expression for M q . In the 
following we will use the Hamiltonian (|43[) in analyzing 
of Anderson localization of the CM of a bright soliton. 
Second order contributions, with respect to the coupling 
to the quasiparticle modes, to the effective Hamiltonian 
(43) are of the order of NVq /n and they can be neglected 
for the parameters of the system used in the present pa- 
per. 



III. QUANTUM DESCRIPTION 

From the point of view of quantum mechanics the clas- 
sical ground state solution §5$ breaks U(l) gauge and 
translation symmetries of the quantum many Hamilto- 
nian [25[ . That is, the quantum Hamiltonian commutes 
with U = e lNe and, in the absence of a disorder po- 
tential, also with the translation operator. In the Bo- 
goliubov description; the 9 and q degrees of freedom ap- 
pear as zero energy modes and, thanks to the Dziarmaga 
approach, we know how to properly describe arbitrarily 
large changes in 9 and q. 

The quantum mechanical version of (|3ip reads 

& = I dzv(z) \Mz-q)\ 2 



IV. ANDERSON LOCALIZATION 

We discuss in more detail the Anderson localization 
in the so called optical speckle potential, as realized e.g. 
in the experiment fl3| . The potential originates from 
the light shifts experienced by the the atoms in the laser 
light detuned from the resonance. In effect, the poten- 
tial V(z) oc a|i?(z)| 2 is proportional to the intensity of 
the local field E(z) and to the atomic polarizability a, 
whose sign depends on the detuning of the external light 
frequency from the atomic resonance. 

Any disordered potential i s completely ch aracterized 
by its correlation functions V(z\) . . . V(z n ) where the 
overbar denotes an ensemble average over disorder re- 
alizations. The average potential value shif ts the ori- 
gin of energy and can always be set to zero, V{z) = 0. 
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Figure 1: Dashed lines: bare potential V(z), solid lines: con- 
voluted potential, i.e. / dz'V(z')\<j>o(z' - z)\ 2 /N. Panel (a) 
for the bare potential amplitude Vo = +0.1 (red detuned laser 
case), panel (b) for Vo = —0.1 (blue detuned laser case). The 
correlation length of the bare potential is ao — 0.28£ where 
£ = 0.02. 
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Figure 2: Panels (a) and (b): eigenstates of the CM of a 
bright soliton; panels (c) and (d): corresponding probability 
density in log scale. The eigenstates correspond to the CM 
momentum P q ~ 10. The red detuned laser case is shown in 
(a) and (c) while the blue detuned one in (b) and (d). The 
inverse localization length is 7 = 23 ± 3 (red detuned case) 
and 7 = 16.5 ± 0.6 (blue detuned case). The parameters of 
the potentials are the same as in Fig. [T] 



The pair correlator can be written as V(z')V(z' + z) = 
VqC(z/(Jo), where Vo measures the potential strength, 
and (To the spatial correlation length. For a gaussian 
random process, higher order correlation functions are 
simple functions of the average and the pair correlator. 
This is no longer the case for non-gaussian potentials that 
require to specify also higher-order correlations. 

An optical speckle potential is a good example of such a 
non-gaussian behaviour. At fixed detuning, the potential 
features either random peaks (the "blue-detuned" case) 
or wells ("red-detuned"). Even after shifting to V(z) = 0, 
the potential distribution is asymmetric (compare Fig.[T]), 
and the importance of odd moments can be probed exper- 
imentally by comparing the blue- and red-detuned cases 
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Figure 3: Panels (a) and (b): eigenstates of the CM of soliton; 
panels (c) and (d): corresponding probability density in log 
scale. The eigenstates correspond to the CM momentum P q sa 
50. The red detuned laser case is shown in (a) and (c) while 
the blue detuned one in (b) and (d). The inverse localization 
length is 7 = 0.27±0.03 (red detuned case) and 7 = 1.8±0.1 
(blue detuned case). The parameters of the potentials are the 
same than in Fig. [1] 



for fixed amplitude |Vo|- The latter is determined by the 
laser strength, and we will use \Vq\ — 8 • 10~ 5 |/i| = 0.1 in 
the following. The bare speckle potential has the pair 
correlation function C(y) = [sin(y)/y] 2 , with a corre- 
lation length that can be as short as 0.28 /im or 
(To = 0.0056 in our units. We shall use this value in the 
following. The CM of the soliton feels, however, not the 
bare potential, but rather its convolution with the soliton 
shape, see Eq. (j4"51) . The convoluted effective potential 
/ dz'V{z')\(j) {z' - z)\ 2 /N (the N factor in the denom- 
inator is due to the normalization of 4>o) is also shown 
in Fig. [TJ While the convolution makes the potential 
smoother it is apparent that it remains quite asymmet- 
ric, thus we may expect that the non-gaussian character 
(in particular non-vanishing odd moments) shows up in 
the properties of the system. For that reason we compare 
the results for both red- and blue- detuned potential of 
similar amplitude. 

The generic properties of Anderson localization in ID 
[3] allow us to expect that all the eigenstates of are 
exponentially localized, i.e., have a typical shape with 
the overall envelope 



|*| 2 ocexp[- 7 (P,)| g -<fo|]. 



(44) 



with qo being the mean position while j(Pq) is naturally 
referred to as the inverse localization length. It depends 
on the eigenenergy of the state E, or writing P q w y/2NE 
on the associated momentum P q . By diagonalizing the 
Hamiltonian (|43[) on a grid, we obtain the wavefunctions 
that are represented in Fig. [5] and Fig. for two signif- 
icantly different energies (momenta). Fig. [5] shows the 
probability densities for the CM of the soliton at rel- 
atively low energies, observe that the exponential enve- 
lope behaviour is visible over several decades. Due to the 
tridiagonal form of the diagonalized matrix on the grid, 
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Figure 4: Inverse localization length as a function of the mo- 
mentum for red-detuned (red-line, solid) and blue-detuned 
(blue-line, dashed) potentials obtained using the transfer ma- 
trix technique. The red dots as well as the blue squares cor- 
respond to the wavefunctions shown in Fig. [2] and Fig. [3] 
Observe the exponential decay of 7 for sufficiently large P q . 

the errors are well under control and the accuracy seems 
not to be limited by double precision arithmetics. Ob- 
serve that the inverse localization lengths obtained for 
red-detuned case and the blue-detuned situation differ 
significantly, stronger localization is observed for the for- 
mer case. 

The situation is quite different at higher energies as 
shown in Fig. [3l Observe that now blue-detuned po- 
tential leads to a much stronger localization. Of course 
the inverse localization lengths at high energies are much 
smaller than those depicted in Fig. fJJ in fact, at suffi- 
ciently high energies 'y(Pq) decays exponentially with P q 
as observed by us before [21J . 

The inverse localization lengths shown as lines in Fig. [4] 
are obtained by a transfer matrix technique [28i ] and quite 
nicely agree with values obtained from exact diagonal- 
ization. Clearly there is a striking difference between the 
two cases of red-detuned and blue-detuned potentials ex- 
emplifying its non-gaussian character and the importance 
of higher moments, in particular third moments. This in 



turn indicates that the application of the celebrated Born 
approximation [H, [2l| which considers the two lowest mo- 
ments only is deemed to fail in our case despite the fact 
that the potential is very weak, smooth and thus, at first 
glance one could naively expect the Born approximation 
to perform quite well. 

With exponentially localized eigenstates, one can now 
consider the dynamics, e.g., the spread of an init ially 
localized wavepacket. As shown by us elsewhere |2l| . 
one can expect an algebraic localization of such a CM 
wavepacket. For realistic parameters, localization occurs 
on a timescale of seconds making the experimental veri- 
fication of the localization feasible. We refer the reader 
to [U for details. 

V. CONCLUSIONS 

Using the Bogoliubov expansion and treating the zero 
modes non-perturbatively, we have shown in detail how 
to obtain the effective quantum Hamiltonian which gov- 
erns the motion of the center of mass of the bright soliton 
in a weak and smooth potential without affecting the soli- 
ton shape. When this potential is of the disorder type one 
may expect to observe Anderson localization of the CM 
motion. The optical speckle potential was considered as 
a realistic example. It turns out that localization proper- 
ties of wavefunctions strongly depend on the sign of the 
potential (red- or blue- detuning). This indicates that, 
even for a weak potential, applicability of the Born ap- 
proximation is limited and the quantitative predictions 
depend on higher correlation functions of the disorder 
potential. Anderson localization of the CM of a bright 
soliton should be experimentally observable. 
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